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Abstract 

The deformation, failure, and low cycle fatigue life of SCS-6/TM5-3 composites are 
predicted using a coupled deformation and damage approach in the context of the 
analytical generalized method of cells (GMC) micromechanics model. The local effects 
of inelastic deformation, fiber breakage, fiber-matrix interfacial debonding, and fatigue 
damage are included as sub-models that operate on the micro scale for the individual 
composite phases. For the laminate analysis, lamination theory is employed as the global 
or structural scale model, while GMC is embedded to operate on the meso scale to 
simulate the behavior of the composite material within each laminate layer. While the 
analysis approach is quite complex and multifaceted it is shown, through comparison 
with experimental data, to be quite accurate and realistic while remaining extremely 
efficient. 


1. Introduction 

Advanced titanium matrix composites (TMCs) continue to tantalize the aerospace research and 
design community. For years these materials have offered seemingly limitless potential for super 
lightweight components with outstanding mechanical properties. However, an array of challenges, such 
as manufacturability, fiber-matrix incompatibility, poor transverse properties, and environmental attack, 
have yet to be overcome and have, in a large part, kept TMCs in the realm of research and on the sideline 
in terms of design. 

An additional challenge related to TMCs involves the development of accurate material and 
structural design and analysis tools. The titanium matrices that are mainly responsible for these 
composites’ attractiveness for aerospace applications undergo inelastic and time/history dependent 
deformation - even during manufacturing. Residual stresses that arise during the cooling of the 
composite from the consolidation temperature, while present in all composites, have a major impact on 
the behavior of TMCs. Debonding at the interface of the fiber and matrix is an additional characteristic 
typically found in TMCs. Further confounding the situation is the fact that most of these phenomena are 
not directly measurable. Clearly, analysis of TMCs requires more than just simple linear elastic 
approaches, such as rule of mixtures equations, which have proven adequate for the ubiquitous polymer 
matrix composites. In order to predict the response of TMCs with any degree of accuracy, more advance 
and complex analysis approaches are needed that account for the characteristics of TMCs described 
above. 
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One micro scale analysis approach that has proven to be quite accurate (in certain situations) for 
TMCs is that of finite element analysis (FEA). Commercial FEA software packages, which allow the 
incorporation of arbitrary inelastic constitutive models for the metal matrix, have been used extensively 
with unit cell approaches to model TMCs (Mall and Nicholas, 1998; Clyne and Withers, 1993). 
However, the accurate inclusion of effects such as residual stresses, fiber-matrix debonding, and fiber 
breakage are difficult. Further, FEA based approaches often require complex boundary conditions to be 
applied to the unit cell. This can make applying different types of loading combinations cumbersome. 
Finally, FEA unit cell approaches require a large number of elements, so if analysis of a number of fibers 
in a certain arrangement is desired, or if the composite being analyzed is part of a larger structural 
problem, the FEA problem can quickly become intractable. 

The approach taken herein employs a micromechanics model known as the generalized method of 
cells (GMC) (Aboudi, 1991, 1995). This model is analytical in nature (as opposed to the numerical FEA) 
and its formulation involves application of several governing conditions in an average sense. This 
averaging renders the model less accurate than FEA at the micro scale due to a decoupling between 
normal and shear field components, but makes it many times more efficient computationally. Unlike 
some analytical models, GMC does provide the local fields in composite materials, allowing 
incorporation of arbitrary inelastic constitutive models for the composite phases as well as other micro 
scale effects such as fatigue damage, fiber breakage, and fiber-matrix debonding. Perhaps most 
importantly, GMC results in closed-form constitutive equations for the composite, which, combined with 
its efficiency, makes the model ideal for inclusion within larger structural analysis approaches to simulate 
the composite behavior at a material point. 

In this investigation, many of the analysis considerations discussed above are exercised. The 
TMC under consideration is silicon carbide/titanium (SiC/Ti). The particular fiber and matrix alloy are 
SCS-6 and Ti-15-3, respectively. This material system has attracted a significant amount of attention 
over the past decade due to its relative ease of manufacture, low cost, and attractiveness as a model 
material system (Arnold and Castelli, 1995) for high-temperature aerospace structural applications (e.g., 
in air breathing engines). GMC has been employed to simulate the tensile deformation, failure, and low 
cycle fatigue (LCF) life of unidirectional SCS-6/Ti-15-3, as well as SCS-6/TM5-3 laminates. The 
inelastic behavior of the Ti-15-3 has been modeled using incremental (time-independent) plasticity 
theory, which is relatively easy to characterize. Local (micro scale) models for longitudinal fiber 
breakage (Curtin, 1991), fiber-matrix debonding (Bednarcyk and Arnold, 2000a), and fatigue damage 
(Arnold and Kruch, 1994) have been included as well. For analysis of the TMC laminate, GMC has been 
embedded within lamination theory, which is used to model the laminate as a structure. In these 
simulations, GMC operates on the meso scale, between the global scale of the laminate and the micro 
scale of the individual phases in the composite. The GMC analysis represents the behavior of the 
composite material at the integration points within each layer of the laminate. The approach is fully 
coupled in that the analysis on each scale affects the other scales as localization and homogenization 
occur repeatedly throughout the application of the simulated applied loading. Further, the fatigue damage 
analysis couples the micro, meso, and global scales as small amounts of damage are applied incrementally 
and after each applied increment, the stress in the laminate is recalculated to allow for the redistribution 
due to the accumulation of damage. Clearly, an equivalent analysis would be (at best) quite cumbersome 
using only FEA. 
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2. Micromechanics model: The Generalized Method of cells 

The micromechanics model employed to simulate the response of the SiC/Ti-15-3 composite is 
the generalized method of cells (GMC) developed by Aboudi (1991, 1995). The geometry of the doubly 
periodic version of GMC is shown in Fig. 1, wherein the microstructure of a periodic material is 
represented by a rectangular repeating unit cell consisting of an arbitrary number of rectangular subcells, 
each of which may be a distinct material. The method assumes a linear displacement field in each subcell 
and imposes continuity of traction and displacement components between subcells in an average sense 
(Aboudi, 1995). This procedure results in a system of linear equations, 

A Ej .-D(e'’ + a,Ar) = Ke 0) 

which are solved for the vector of all subcell strains, E s , to form strain concentration equations, 

s, = Ac + D(ef +a,AF) (2) 

In the above equations, A , D , and K are matrices containing information on the repeating unit cell 
geometry and the properties of the constituent materials, E is the vector of the global (unit cell) strain 

components, ef is the vector of all subcell inelastic strain components, a s is the vector of all subcell 

coefficients of thermal expansion (CTEs), and A and D are strain concentration matrices. From eq. (2) 
the effective thermo-elastoplastic constitutive equations for an arbitrary periodic material, that are 
required to simulate the response of the composite material to external loading, can easily be determined 
(Aboudi, 1995). Clearly, in this original formulation of GMC, the subcell strains, £ 5 , serve as the basic 
unknown quantities. Further, in solving eq. (1) for these quantities, a great deal of computational effort 
may be required as the matrix A becomes large as the number of subcells in the analyzed unit cell 
increases. 

The recent reformulation of the GMC equations (Pindera and Bednarcyk, 1999; Bednarcyk and 
Pindera, 2000) employs subcell stresses (rather than strains) as basic unknowns and develops mixed 
concentration equations for the heterogeneous material. Due to GMC’s inherent lack of normal-shear 
field coupling, significantly fewer unique subcell stress components are present compared to subcell strain 
components (see Pindera and Bednarcyk (1999) and Bednarcyk and Pindera (2000) for details). Tracking 
only the unique subcell stress components gives rise to, 

GT=r-f'-f p , (3) 

which replaces eq. (1). Here, G is a matrix containing information on the repeating unit cell geometry 
and the properties of the constituent materials, T is the vector of all unique subcell stress components and 
f , f ' , and i p are vectors containing global total strains, subcell thermal strains, and subcell inelastic 
strains, respectively. Solving eq. (3) provides mixed concentration equations, 

T = Gs + r AT + 0 , (4) 

where the thermal and inelastic terms are included in T and O , respectively. These are mixed 
concentration equations because they provide local stresses in terms of global strains; G (which is the 
inverse of G ) is the mixed concentration matrix. Clearly, since there are far fewer unique subcell stress 
components than strain components, a significant improvement in efficiency is gained by employing the 
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reformulation of GMC. Thus, by exploiting GMC’s lack of shear coupling, the reformulation represents 
an ultra-efficient multi-axial micromechanics analysis tool. Direct comparison of execution times 
resulting from the reformulated and original versions of GMC indicates a striking speed-up associated 
with the reformulation, particularly as the number of subcells becomes large. For example, for a 12x12 
subcell repeating unit cell, the reformulated version of GMC has been shown to be more than 19,000 
times faster (Bednarcyk and Arnold, 2000a). As demonstrated by Pindera and Bednarcyk (1999), the 
original formulation and reformulation of GMC yield identical results, thus the speed-up comes with no 
loss of accuracy. Clearly, utilization of the reformulation is crucial to the realization of reasonable 
execution times for application of GMC to detailed composite microstructures and when embedded 
within lamination theory to represent the material response at an integration point. In fact, most of the 
results generated for this study could not have been generated (due to exorbitant execution times) prior to 
the reformulation of GMC. 

In order to simulate the response of composite laminates, the reformulated version of GMC has 
been embedded within classical lamination theory, as shown in Fig. 2. The basic thermo-elastic 
lamination theory equations are operable (see Jones, 1968) with the additional terms that represent the 
effects of inelastic strains. The overall laminate constitutive equation can then be written, 


"a b" 

V 


"N" 


~N r " 


" N"' ~ 



= 


+ 


+ 


B D 

K 


M 


_M r _ 


_M"' 


where A , B , and D are the extension, coupling, and bending stiffness matrices, £° is the vector of 
midplane strain components, K is the vector of midplane curvatures, N and M are the mechanical force 

and moment resultant vectors, N r and M 7 are the thermal force and moment vectors, and N w and M"' 
are the inelastic force and moment vectors (for additional details see Bednarcyk and Pindera (1996) and 
Herakovich ( 1 998)). 

As indicated in Fig. 2, while the overall laminate response is dictated by eq. (5), the response of 
the composite material at each Gauss integration point within each layer is determined via GMC. The 
effective stiffness and CTEs of each layer, once properly transformed to the global (laminate) coordinate 
system via standard rotation equations, allow the determination of the A, B, and D matrices and the 
thermal force and moment vectors (for a given temperature change). In the absence of inelasticity, the 
formulation would then be complete. Given laminate scale applied loading in terms of an appropriate 
combination of mechanical force and moment resultants and midplane strains and curvatures along with a 
temperature change, the global response of the laminate can be determined using eq. (5), However, in the 
presence of inelasticity, the model must constantly localize and homogenize as the loading is applied in 
an incremental fashion. That is, the laminate scale loading results in a stress and strain field at the 
integration points within each layer. These fields are transformed (via rotation) and passed to GMC, 
which further localizes the fields from the level of the unit cell to each individual subcell. Then, on this 
micro scale, the inelastic strains are determined via an arbitrary inelastic constitutive model, and 
homogenized to yield the cell (meso scale) inelastic strains. This inelastic strain field, once transformed 
to the laminate coordinate system, is representative of a particular integration point in a particular layer of 
the laminate. Finally, the inelastic field is integrated through the thickness of each layer and the inelastic 
force and moment vectors are determined via the appropriate lamination theory equations. As stated, this 
localization/homogenization process must occur at each increment of the applied loading. Further, for a 
local inelastic constitutive model that includes a consistency condition (such as the incremental plasticity 
theory employed herein), a number of equilibrium iterations are required at each loading increment to 
ensure that the consistency condition is met. 
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3. Constituent response 


The composite material examined in this study consists of a Ti-15-3 matrix reinforced by 
continuous SCS-6 silicon carbide fibers. In order to determine the deformation response of the SiC Ti- 
1 5-3 composite, it is first necessary to incorporate sub-models for the behavior of the individual phases. 
Towards this end, the SCS-6 fiber was treated as an elastic, isotropic material with the temperature- 
dependent properties given in Table 1. The elastoplastic constitutive response of the Ti-15-3 matrix was 
modeled using Mendelson’s (1968) time-independent incremental plasticity theory. This theory pro\ides 
the inelastic strain increment for an isotropic material (with some manipulation of Mendelson’s (1968) 
equations). 


de™ = 




3Ge. 


et J 


( 6 ) 


where G is the elastic shear modulus, e t , is the equivalent modified total strain deviator, and <?' are the 
components of the modified total strain deviator (see Mendelson (1968) for details). The final term 
appearing in eq. (6), rr^, , is the effective stress. This term is determined from the post-yield hardening 

response of the material. Typically, the post-yield response of the material is taken to be bilinear (see 
Bednarcyk and Pindera (2000), Williams and Pindera (1997)), in which case the effective stress is given 


by, 


cr eff =c7,.+- 


EH 


E-H 


'efi 


(7) 


where is the effective plastic strain and H is the secondary, or post-yield, slope of the material’s 

uniaxial stress-strain curve (approximated as bilinear). The present implementation of the incremental 
plasticity theory has been extended to admit piecewise linear post-yield behavior. This allows a more 
accurate representation of the actual inelastic material behavior as the response can be broken into an 
arbitrary number of linear sections. Then, in the / th linear section, the effective stress is given by, 

< 8 > 


where Y t and £, p are points of stress and plastic strain along the piecewise linear representation of the 
uniaxial stress-strain curve, and H\ p is the local slope of the curve. The plastic strain points of the 
piecewise linear representation are calculated from total strain points (which are actually input to the 
analysis code) from, 


£ p 



(9) 


where £ i are the input strain points. 

The temperature-dependent elastic properties for the Ti-15-3 matrix are given in Table 2. The 
inelastic properties, in terms of yield stress and additional stress-strain point pairs, are given in Table 3. 
Figure 3 shows the characterization of the Mendelson incremental plasticity model for the Ti-15-3 matrix 
material (Tuttle and Rogacki (1991)). Clearly, the point-wise nature of the implementation has enabled 
the model to reproduce the temperature-dependent uniaxial stress-strain response of the material quite 
well. Linear interpolation is used to allow the properties to span the entire temperature range indicated. 
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4. Unidirection al Composite Tensile Response 

In order to model the deformation, failure, and fatigue lives of SiC/Ti-15-3 laminates, two 
additional local effects, which pertain to the response of unidirectional SiC/Ti-15-3, must be considered. 
The first is interfacial debonding between the fiber and the matrix due to a local tensile normal stress. 
The second is breakage of individual fibers due to a local axial stress. 

It is well known that unidirectional SiC/Ti composites exhibit fiber-matrix debonding when 
subjected to off-axis loading (see Johnson et al. (1990), Mall and Nicholas (1998), and Bednarcyk and 
Arnold (2000b)). In order to capture even the qualitative nature of the off-axis SiC/Ti-15-3 tensile 
response, a simulation must account for this interfacial debonding (which typically causes a noticeable 
“knee” in the stress-strain curve). A model for interfacial debonding known as the evolving compliant 
interface (ECI) model was recently developed and implemented within GMC by Bednarcyk and Arnold 
(2000b). The ECI model employs the concept of a flexible interface (Jones and Whittier, 1967), wherein 
a discontinuity in the normal or tangential displacement component at an interface, /, is permitted. These 
discontinuities are taken to be proportional to the appropriate stress component at the interface such that, 


[u„] = R„<j, 
[u,] =11,(7, 


( 10 ) 


where [//„ ] and [ft,] are the normal and tangential displacement discontinuities at the interface, cr n | / 

and a )\ f are the normal and tangential stress components at the interface, R n and R t are empirical 

debonding parameters that represent the effective compliance of the interface, and cr n DB and <t' dh are the 
normal and tangential debond stresses (i.e. strengths) of the interface. After debonding occurs, the ECI 
model allows interfacial debonding to progress through unloading of the stress at the interface by treating 
the debonding parameters as time-dependent. The form of this time-dependence is, 


*(') = A 


exp 


t-t 


DB 


B 


-1 


t>t 


DB 


( 11 ) 


where A and B are empirical constants specific to the interface, t is the time, and t DB is the time at 
which debonding occurs, 

Bednarcyk and Arnold (2000b) used GMC, with the ECI model functioning as a sub-model, to 
simulate the transverse tensile and creep response of SCS-6/TIMETAL 21 S composites. In the present 
study, the ECI model (still functioning within GMC) was re-characterized for application to SCS-6/Ti-15- 

3. Further, Bednarcyk and Arnold (2000b) employed a 26x26 subcell repeating unit cell to represent a 
single refmed fiber and allowed debonding to progress around the fiber during simulated loading. Herein 
the employed repeating unit cell contains 28 individual fibers, each with its own interfacial bond strengths 
and A parameters. The bond strengths and parameters employed in the present study are given in Table 

4. The parameters in Table 4 associated with normal debonding were characterized based on a 
comparison of GMC predictions with experiment for a room-temperature transverse tensile test on a 35% 
fiber volume fraction SCS-6/Ti-15-3. This characterization is shown in Fig. 4. Note that residual stresses 
from fabrication were included in the model simulations via application of a globally stress-free cool 
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down (from an assumed stress-free temperature of 600 °C)' prior to application of the simulated tensile 
loading. The composite’s simulated failure was accomplished using a global transverse maximum stress 
failure criterion of 429 MPa. 

Due to the unavailability of experimental shear test results for SCS-6/Ti-15-3 and the need for 
capturing multiaxial effects, the parameters in Table 4 associated with tangential debonding weie 
assembled using the normal debonding parameters as a starting point. Alterations were then made to 
better capture the qualitative nature of the [±45]> s laminate deformation response. This correlation is 
shown in the next section. 

The breakage of individual fibers that leads to the overall failure of longitudinally loaded 
unidirectional SiC/Ti composites also must be addressed via a sub-model. It is well known that the rule 
of mixtures type failure models, which assume that the composite’s fibers fail when the average fiber 
stress reaches the average fiber strength, tend to over predict the strength of SiC/Ti composites 
significantly (e.g., Brindley and Draper (1993). In reality, early breakage of individual lower strength 
fibers decreases the global stiffness of the composite and redistributes additional loads to the intact fibers, 
resulting in lower axial strengths. Curtin’s (1991) effective fiber breakage model was employed to 
account for these effects. The Curtin model was recently incorporated within GMC by Bednarcyk and 
Arnold (2000a) and applied to model the longitudinal tensile response of SCS-6/TIMETAL 21 S 
composites. Using fiber strength statistics combined with a shear-lag analysis, Curtin (1991) developed a 
general equation describing the stiffness degradation of an effective fiber that represents all fibers in the 
composite as the composite is loaded longitudinally. The effective fiber elastic modulus is given by, 



1 + exp 


f £ mech V " 4 1 


( 12 ) 


where E f is the original fiber modulus, £™ ch is the fiber longitudinal mechanical strain, m is the fiber 

Weibull modulus, and a c = (2<t” t L 0 jdj ' ^ . The term cr 0 is the mean fiber strength, r is the 

frictional sliding resistance between the fiber and the matrix, L 0 is the fiber gauge length, and d is the 
fiber diameter. The model predicts composite failure when the stress in the effective fiber reaches a 
maximum, 

< 7 ™ x = j<t c x l/(m+1) [l + exp(-x)] (13) 

where x is the first positive integer that satisfies, 

l + [l-(»7 + l)jt']exp(-x) = 0 (14) 

One limitation of the Curtin model (revealed by the off-axis laminate results shown later) is its 
uniaxial nature. However, for unidirectional SiC/Ti composites loaded longitudinally, it is still quite 
capable of predicting stiffness degradation and failure when employed within GMC (Bednarcyk and 
Arnold, 2000a). 

Figure 5 compares the longitudinal tensile deformation and failure of SCS-6/Ti-15-3 composites 
as predicted by GMC with experimental results. The employed Curtin model parameters are: 


1 The 600 °C stress-free temperature was selected based on correlation with experiment for the room temperature 
41% SCS-6/Ti- 1 5-3 tensile results shown in Fig. 5. 
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d = 142 pm, L 0 = 12.7 mm, m - 17.0, cr 0 =3.85 GPa (from Draper et al., 1992), and r 0 = 14 MPa 

(from Bednarcyk and Arnold, 2000a). Since the implementation of the Curtin model did not include 
temperature-dependent model parameters, in order to account for the fiber strength degradation due to 
elevated temperature, for the 427 °C predictions, the <7 0 value was reduced by 3.5% (Bednarcyk and 
Arnold, 2000a). Since the Curtin model represents the behavior of all fibers within the composite as one 
effective fiber, a simple 2x2 GMC repeating cell containing a single fiber was employed. 

The room temperature 41% SCS-6/Ti-15-3 results shown in Fig. 5 were employed to characterize 
the composite’s stress-free temperature (600 °C). For this stress-free temperature, the deformation 
predicted by the model agrees well with experiment. The 427 °C results represent pure predictions. 
Examining the 427 °C predictions for fiber volume fractions of 41%, 35%, and 15%, it is clear that the 
model under predicts the deformation somewhat, exhibiting a greater amount of inelastic deformation 
compared to experiment. This indicates that the residual stresses, as simulated by the model using the 
stress-free temperature of 600 °C, are inexact (Goldberg and Arnold, 1999). This is expected as the 
matrix constitutive model is time-independent and thus does not allow any relaxation of residual stresses. 
The model still does a reasonably good job of predicting the composite UTS for these three fiber volume 
fractions. The results for the fiber volume fraction of 25%, on the other hand, show a large discrepancy 
between model and experiment. The cause of this discrepancy is not known at this time, however, the 
experimental results exhibit much better agreement with model predictions for a composite with a fiber 
volume fraction of 28%. “ 

5, Laminate Tensile Response 

Using the information garnered from the tensile response of the unidirectional SCS-6/Ti-15-3 
composites, the lamination theory implementation of GMC can be used to predict the tensile response of 
SCS-6/Ti-15-3 laminates. For the laminate, the coupled analysis spans three levels of scale. The global 
or structural scale problem is addressed using the lamination theoiy equations. The meso scale problem is 
represented by GMC, analyzing the behavior of the composite material at the integration points within 
each layer of the laminate. Finally, the micro scale consists of the individual fiber and matrix subcells 
within the GMC repeating unit cells at each integration point. On the micro scale, the equations of 
incremental plasticity, as well as the ECI debonding and Curtin fiber breakage models, are operative. 
This multi-scale analysis is coupled in the sense that during each increment of the applied loading and 
each equilibrium iteration, the entire scale is spanned through the localization and homogenization 
provided by the models. Further, the predictions of each model on each scale affect the predictions of the 
models functioning on the other scales. 

The GMC repeating unit cells employed for each integration point in each layer of the laminate 
were identical to that employed for the transverse composite simulations. This unit cell consists of 28 
fibers, which were assigned debonding parameters according to Table 4. In addition, the Curtin model 
was used to simulate the longitudinal stiffness degradation and breakage of the fibers. Since the Curtin 
model simulates many fibers as one effective fiber, identical Curtin model parameters (given above in 
Section 4) were employed. 

Figure 6 compares model results with experiment (Lerch and Saltsman, 1993) for the deformation 
and failure of [0/90] 2s , [±45] 2s , and [±30] 2s laminates along with the longitudinal and transverse (i.e., [0] s 
and [90] 8 ) results shown earlier. As mentioned previously, it was necessary to utilize the [±45] 2s results to 
determine the shear debonding parameters for the ECI model (see Table 4). Consequently, the 
deformation results of only the [0/90] 2s and [±30] 2s laminate should be viewed as predictions. Further, a 
repeating unit cell level maximum shear stress failure criterion of 220 MPa (chosen based on correlation) 
was employed in the simulations. Clearly, the model agrees quite well with the laminate experiments. 
Figure 7, however, shows that the GMC model prediction is somewhat quantitatively inaccurate for a 
[±60] 2s laminate. While the experimental stress-strain curve for this laminate falls below the curve for the 
[90] g composites, GMC predicts (as one might expect) that the [±60] s laminate curve falls above that of 
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the [90] 8 composite. Qualitatively, the predicted curve is quite good as it resembles the experimental 
curve closely. The specific stages of the laminate deformation, involving strong bonding with elastic 
behavior, followed by debonding, and finally by large-scale inelastic behavior, are well represented (see 
Bednarcyk and Arnold, 2000b). The quantitative discrepancy is likely due to inadequacies of the GMC 
model in terms of shear (i.e., its lack of normal-shear field coupling). Further, the shear debonding 
parameters could be altered to provide better correlation with experiment in Fig. 7, however, this would 
affect the agreement between model and experiment for the [±45] 2s laminate in Fig. 6. 

6. Laminate Low Cycle Fatigue Response 

With the illustrated ability of GMC to model the deformation and static failure response of 
SiC/Ti-1 5-3 laminates, the LCF life of the composite can now be addressed using an appropriate local 
fatigue damage model. Towards this end, the multiaxial, isothermal, continuum damage mechanics 
model of Arnold and Kruch (1994) has been employed. This model was previously incorporated within 
GMC by Wilt et al. (1997). It assumes a single scalar internal damage variable, D, that has a value of 
zero for undamaged material and one for a completely damaged (failed) material. A summary of the 
damage model equations is given in the appendix. 

The implementation of the damage model within GMC has been performed on the local scale, 
thus damage evolves in a given subcell based on the local stress state and number of cycles. For a given 
damage level, the stiffness of the subcell is degraded by (1 - D), while the inelastic strains are calculated 
using an effective stress that is increased by the factor 1/(1 - D ). Further, the GMC implementation 
allows the application of a local damage increment, A D , and then calculates the number of cycles, N, 
required to achieve this local increment of damage. This approach allows the model to determine the 
stress state in the composite, identify the subcell that will reach the desired damage level in the fewest 
cycles, apply that number of cycles, and calculate the damage that arises throughout the remainder of the 
composite. Then the composite can be reanalyzed and a new stress state determined based on the new 
damage level throughout the composite. In this way the local and global stress and damage analyses are 
coupled. As the damage in the composite evolves, the stress field in the composite is redistributed, which 
then affects the evolution of damage. 

The damage model was characterized for the isotropic Ti-15-3 matrix at 427 °C. Given the fact 
that the specific material is isotropic, the parameters co u , cOj , , Q ) m , t ] u , , and i] m are set equal to one, 

and the model summarized in the appendix reduces to the well-known ONERA NLCDM model 
(Chaboche and Lesne (1988), Lemaitre and Chaboche (1990)). The parameters /?, a, (T ^ . A/ 0 ,and <7„ 
were selected based on correlation with matrix data for fully reversed stress controlled fatigue tests (i.e., 
the mean stress, a = 0 and R = (r min /<T mas - -1 ). This correlation is shown in Fig. 8. The parameters 

b and b' determine the mean stress dependence of the damage model parameters M and a n such that, 

M = M 0 (\-ba) (15) 

The correlation for selecting b and b' , based on experiments with a = 0 , is shown in Fig. 9. Note that 
the model results in Figs. 8 and 9 are purely from the damage model, GMC was not involved in these 
calculations. 

In the case of the elastic fiber, a simple strength degradation damage model (which is intended to 
account for interfacial and fiber damage due to repeated cyclic loading), in accordance with the study by 
Wilt et al. (1997), was implemented for the SCS-6 fiber. This type of model was introduced rather than 
using the previous fatigue damage model since no marked change in the stiffness of the fiber due to cyclic 
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loading is experimentally observed. This implementation involved degradation of the Curtin model fiber 
strength parameter, <r 0 , with the number of applied fatigue cycles, N. Correlation, performed using 

longitudinal LCF life results, yielded the following degradation equation for <x 0 , 

a 0 = 3.72 GPa N< 10,000 

< 7 0 = 11.8 — 2.02 log(A r ) GPa 1 0,000 < TV <23,263 (17) 

cr 0 = 2.98 GPa N> 23,263 

This form and these values were chosen based on correlation using the [0] 8 experimental data shown in 
Fig. 10. 

Using the above-described damage models for the fiber and matrix, the Mendelson incremental 
plasticity for the matrix inelastic behavior, the ECI model for the normal and tangential debonding at the 
fiber-matrix interface, and the Curtin model for the longitudinal breakage of the fiber, GMC was 
employed to predict the LCF life of [±30]2 S and [±45]2 S laminates, as well as the transverse [90] 8 
composite at 427 °C. The results are shown in Fig. 10. Generation of the shown model predictions was 
quite intense from a computational standpoint as each of the two integration points in each of the three 
layers of the laminate was represented by a GMC unit cell consisting of 28 fibers (112 subcells). A 
globally stress-free cool down from 600 °C to room temperature, followed by a globally stress-free heat 
up to 427 °C, was applied prior to application of the simulated applied loading cycle. Recall that due to 
the presence of plasticity, the loading was applied incrementally with equilibrium iterations required for 
each increment. The loading cycle was then repeated for each applied local damage increment until 
complete failure of the laminate, resulting from an accumulation of local and meso scale failures, 
occurred. Typically, a local damage increment of 0.2 was employed, the loading cycle was broken into 
240 total increments, an average of 20 equilibrium iterations were required for convergence, and 
execution of the model (for each <J m value) took in the range of 10 to 20 minutes 2 on a 700 MHz Pentium 

III personal computer. Considering the complexities of this multiscale approach and the number of 
variables the simulations required, the authors regard the model’s execution as being quite efficient, 
especially when compared to the alternatives available for such an analysis. 

It is clear from Fig. 10 that GMC, with its incorporated micro scale models, does quite a good job 
of matching with the experimental data (Lerch (1990) and Lerch et al. (1993)) for this complex problem. 
The [90]g, [±45] 2s , and [±30] 2s model results are pure predictions, while, as described above, the [0] 8 
results have been correlated below a stress level of 1.04 GPa. The agreement with experiment for the 
[±45] 2s laminate is remarkable. For the [±30] 2s laminate, the agreement appears to be good in the higher 
life (low stress) range. However, as the predicted cycles to failure decrease, the stress increases to a much 
greater degree than in the experimental data. That is, the typical “step” in the predicted S-N curve is too 
severe, and the life is over predicted for the higher stresses. This discrepancy is likely due (in some 
degree) to the uniaxial nature of the Curtin fiber breakage model since the Curtin model is dependent only 
upon the axial stress in the effective fiber, which is decreased significantly in the [±30] 2s laminate 
compare to the [0] 8 composite. Whereas, in reality, the off-axis stresses, which are significant in the 
[±30 ] 2s laminate, should contribute to the fiber breakage. In the [±30] 2s laminate tensile response, this 
problem was overcome by employing a meso scale maximum shear stress criterion (see Fig. 6). 
However, this criterion activates only for the highest applied stress, leading to a life of only 1 cycle in Fig. 
10. This is clearly not a problem for the [±45] 2s laminate, whose behavior is more matrix and interface 
dominated. Additional potential explanations for the discrepancy in the [±30] 2s results include increased 
fiber damage due to the off-axis fiber orientation and inaccuracy of the matrix stress field due GMC's 
lack of normal-shear coupling. 


2 Note that the executed cases were not consciously optimized for speed. 
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In contrast to the [±30] 2s laminate prediction is the [90] g composite prediction, which appears to 
be in good agreement with experiment for the lower life range, but poor agreement in the higher life 
range. The predicted S-N curve is clearly too flat. The [90] 8 composite is dominated to the greatest 
extent by the matrix and interface behavior, and thus this discrepancy is likely due to inaccuracies 
associated with the matrix constitutive model (time-independent plasticity) and the ECI debonding model. 
In addition, once complete debonding of the fiber-matrix interfaces in the composite has occurred, it is as 
if the composite contains holes (rather than fibers), which give rise to stress concentrations in the matrix 
and lower fatigue life. Such stress concentrations are not well represented by GMC due to the averaging 
performed in the method's formulation. This would tend to cause GMC to over predict the life of the 
[90] 8 composite as illustrated in Fig. 10. Considering the complexities associated with predicting the LCF 
life of these laminates along with the approximate or effective nature of many of the modeling approaches 
employed, the overall agreement between model and experiment exhibited in Fig. 10 must be considered 
to be very good. 

7. Conclusion 

A fully coupled deformation and damage approach to modeling the response of composite 
materials and composite laminates has been presented. It is based on the analytical GMC 
micromechanics model, which provides closed-form constitutive equations for composite materials as 
well as the micro scale stress and strain fields in the composite phases. The provided constitutive 
equations allow GMC to function within a higher scale structural analysis to represent a composite 
material point, while the availability of the micro fields allows the incorporation of lower scale sub- 
models to represent local phenomena in the fiber and matrix. Further, GMC’s formulation performs 
averaging when applying certain governing equations, such that some degree of micro scale field 
accuracy is surrendered in favor of extreme computational efficiency, rendering the method quite 
attractive as the centerpiece in a multi scale structural analysis. 

This approach has been applied to simulate the deformation, static failure, and LCF life of 
SiC/Ti-1 5-3 composites. The important effects of residual stresses, matrix inelasticity, fiber breakage, 
and fiber-matrix debonding have been incorporated within GMC through appropriate sub-models and 
load histories. For the LCF calculations, an appropriate continuum damage fatigue model was 
incorporated on the micro scale as well. For simulations of composite laminates, GMC, with all its sub- 
models intact, was embedded within lamination theory, which served as the global or structural scale 
model. The approach then spanned three fully coupled levels of scale for each iteration of each increment 
of the applied simulated loading. Results from the GMC model were compared with experiment and 
exhibited good agreement for the defonnation, static failure, and LCF life of the SCS-6/Ti-15-3 
composite materials and laminates. This type of multiscale analysis of composites and composite 
structures shows promise for future design and optimization tools. 
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Appendix: Fatigue Damage Model 


For additional details on the fatigue damage model, see Arnold and Kruch (1994) and Wilt et al. 
(1997). The evolution of damage is governed by, 




1 - D 


dN 


(Al) 


where N is the number of cycles at the current stress state ( cr k ), D k and D k X are the amount of damage 
after the current and previous damage increments, respectively, and a is a function of the current stress 
state, 

K) 


i G n 
a - 1 —a — 


. (<•>.) 


(A2) 


and ( ) are Macauley brackets. cr /; and o u are the uniaxial fatigue limit and ultimate strength (see eq 
16 and Table 5), respectively, and O y , and O,, are the fatigue limit and static fracture surfaces, 
respectively, 

1 max max 


O fl = 

* 2 /, 


0 = 1 - 



(')- CT s('«)]-l 

(A3) 

max _ 

[^(O] 

(A4) 


and the normalized stress amplitude is defined as, 

1 max max 


F_. =■ 


ax max „ r , , \-i 

* 


(A5) 


When (O„) = 0, static fracture (complete local failure) is indicated, thus a cannot be 

undefined. ^O^ = 0 indicates that the current stress state is below the fatigue limit and a is set to one. 

This then represents a special case when integrating eq. (Al) that will be considered separately. The t 0 
and t terms in eqs. (A3 - A5) are the time at the beginning of the current load cycle and some time 

max 

during the current load cycle, respectively. The operator 


indicates that the maximum value of the 


expression to the right up to time / of the current load cycle should be taken. 

The general form for the terms F, ^ , F [ n t , and can be expressed as, 


F u = 


V< i 


H.-'H 
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(A 6) 
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where, 


I 1 =d,d l S Jl S t ,-{d l <l J S, J ) 1 (AT) 


The current deviatoric stress components are S tj = <J tJ ~~<y mm S IJ , and r/ ( are the components of the 

vector defining the preferred direction in a transversely isotropic material. co { } and // ( ( represent the 

ratios of longitudinal to transverse normal and shear stresses, respectively, for a transversely isotropic 
material and are equal to one for the current isotropic material (Ti-15-3) application. 

For the case in which the current local stress state, cr k , is above the initial fatigue limit, the 

number of cycles, N, required to damage a subcell from a level of D k] to D k is obtained by integrating 
eq. (Al) and is given by, 


N 


{[■ - (i - d, y-' ]' " - [1 - (i - d,., y J" } 


(A8) 


When the local stress level is below the fatigue limit, a = 1 , integration of eq. (Al) yields, 


N=- 


log 


.-(i-D,r]-iog[i-(.-A.r] 

FJUfi + l) 


(A9) 


Equations (A8) and (A9) can also be solved for D k in order to determine the current amount of damage 
developed for a given number of cycles and previous state of damage, D k x . Further, to determine the 
remaining cycles to failure for a given previous state of damage, D k can be set equal to one in eqs. (A8) 
and (A9). 
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Table 1 . Temperature-dependent elastic properties for the SCS-6 fiber (Bednarcyk and Arnold, 2000a). 


T(°C) 

E (GPa) 

V 

a (xlO' 6 /°C) 

21. 

393. 

0.25 

3.56 

316. 

382. 

0.25 

3.73 

427. 

378. 

0.25 

3.91 

538. 

374. 

L_ 0.25 

4.07 

860. 

368. 

0.25 

4.57 


Table 2. Temperature-dependent elastic properties employed for the Ti-1 5-3 matrix (Rosenberg, 1 983). 


T(°C) 

E (GPa) 

V 

a (xlO' 6 /°C) 

23.9 

89.52 

0.3 

8.33 

204.4 

86.38 

0.3 

8.70 

315.6 

78.89 

0.3 

9.13 

482.2 

74.72 

0.3 

9.63 

565.6 

71.12 

0.3 

9.88 

648.9 

65.70 

0.3 

10.23 


Table 3. Piecewise linear, temperature-dependent inelastic data for the Ti-1 5-3 matrix determined via 
correlation with experiment (see Fig. 3). 


T = 23.9 °C 

T = 204.4 °C 

T = 315.6 °C 

T = 482.2 °C 

T = 565.6 °C 

T = 648.9 °C 

<7 V = 710.2 MPa 

<7 y - 566.8 MPa 

<J V = 509.5 MPa 

<T r = 413.0 MPa 

< 7 v = 208.2 MPa 

<j y = 83.43 MPa 

Y 

(MPa) 

S, (%) 

Y 

(MPa) 

e, (%) 

Y 

(MPa) 

£, (%) 

Y 

(MPa) 

£, (%) 

Y 

(MPa) 

E x (%) 

Y 

(MPa) 

S, (%) 

786.0 

0.963 

635.0 

0.823 

572.3 

0.799 

499.9 

0.778 

239.9 

0.379 

142.0 

0.303 

799.8 

1.09 

663.3 

0.999 

609.5. 

1.00 

517.8 

0.995 

252.4 

0.505 

155.1 

0.480 

813.6 

1.51 

696.4 

1.50 

641.9 

1.48 

518.5 

1.50 

255.1 

1.01 

163.4 

0.976 

820.5 

2.01 

717.1 

2.18 

662.6 

2.02 

518.8 

2.00 

259.3 

2.04 

165.5 

1.97 
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Table 4. ECI model parameters of the 28 fibers in the composite GMC repeating unit cell. 


Fiber Number 

<y'm (MPa) 

A„ (GPa' 1 ) 

<r' m (MPa) 

A, (GPa') 

i 

6.9 

14.5 

6.9 

14.5 

2 

20.7 

13.1 


13.1 

3 

41.4 

11.6 

34.5 

11.6 

4 

69.0 

10.2 

48.3 

7.25 

5 

82.7 

2.90 

69.0 

1.45 

6 

138. 

0.725 

148. 

1.02 

7 

159. 

0.435 

159. 

1.02 

8 

186. 


169. 

0.870 

9 

207. 


179. 

0.870 

mm 

262. 



0.725 

ii 

269. 



0.725 

12 

276. 

0.218 

210. 

0.580 

13 

290. 

0.174 

221. 

0.580 

14 

303. 

0.174 

231. 

0.435 

15 

324. 

bbksbbhi 

241. 


16 

338. 


252. 

0.145 

17 

359. 

0.174 

262. 

0.145 

18 

386 


272. 

0.116 

19 

413. 


283. 

0.101 

20 

441. 

0.131 

293. 

0.0870 

21 

469. 

0.0870 

303. 

0.0580 

22 

496. 

0.0870 

314. 

0.0435 

23 

552. 


324. 

0.0290 

24 


0.0435 

Strong 

— 

25 

p— 

0.0435 

Strong 

— 

26 

758. 

0.0261 

Strong 

- 

27 

1000. 

0.0131 

Strong 

— 

28 

Strong 


Strong 

- 
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Table 5. Damage model parameters for the isotropic Ti-1 5-3 matrix. 


parameter 

value 

Units 

b 

0.0 


b' 


- 


1.0 

- 

(Oji 

1.0 

— 


1.0 

— 

n„ 

1.0 

- 


1.0 

— 

>L 

1.0 

— 

P 

1.0 

- 

a 

0.0022 

- 

a \ 

207. 

MPa 

M o 

12.4 

GPa 


725. 

MPa 
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Doubly Periodic GMC Analysis 


Figure 2. Schematic showing the lamination theory geometry with GMC embedded to represent the 
behavior of the composite material at the through-thickness integration points. 
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a Experiment 
— Plasticity Model 


23.9 °C 


204.4 °C 


315.6 °C 


482.2 °C 


kAAAAAAAAAAAAAAAAAAAAAAAAAAA 


565.5 °C 



Figure 3. Charactenzation of the Mendelson incremental plasticity model for the Ti-15-3 matrix at 
various temperatures. These data represent an applied strain rate of 10‘ 4 /s (Tuttle and Rogacki, 
1991). 
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Strain 


Figure 4. Transverse tensile response of 35% SCS-6/Ti-15-3 at room temperature - characterization of 
the ECI model parameters. These data represent an applied strain rate of 10' 4 /s. The curves’ 
end points indicate failure. 
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2000 



A ° Vf = 41% 
Vf = 35% 

Vf = 25% 
Vf = 15% 


o Experiment - 23 C 
GMC Model - 23 C 
a Experiment - 427 C 
— GMC Model- 427 C 


0.000 0.001 0.002 0.003 0.004 0.005 0.006 0.007 0.008 0.009 0.010 

Strain 


Figure 5. Comparison of GMC model results (employing the Curtin fiber breakage model) with 
experiment for the longitudinal tensile response of SCS-6/Ti-15-3. The 23 °C (Vf = 41%) 
results were employed to select the composite’s stress-free temperature of 600 °C. The 
curves’ end points indicate failure. 


NASA/TM— 2001-211343 


24 





Figure 6. Comparison of GMC model results and experiment for the tensile response of SCS-6/TM5-3 
laminates at 23 °C. The fiber volume fraction for the [0] 8 composite was 41%, for all others 
it was 35%. The deformation of the [0/90] 2s and [±30] 2s laminates are pure predictions while 
the remaining model results included some correlation. The curves’ end points represent 
failure in all cases except the [±45] 2s laminate, which experienced a run-out. 
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Strain 

Figure 7. Comparison of GMC model results and experiment for the tensile response of 35% SCS-6/Ti- 
15-3 laminates at 23 °C. 
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Maximum Stress, G m (MPa) 



Figure 8. Correlation for the selection of damage model parameters for the Ti-15-3 matrix at 427 °C. 
These data represent fully reversed tests (cr = 0, R = — l) . 
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Figure 9. Correlation for the selection of damage model parameters for the Ti-15-3 matrix at 427 °C. 
These data represent applied tensile stress states (ct ^ 0, R = 0) . 
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[+ 45 /- 45] 2 


1,000 10,000 

Cycles to Failure, N f 


100,000 


1 , 000,000 


Figure 10. Comparison of model LCF life predictions (S-N curves) with experiment for 35% SCS-6/Ti- 
15-3 at 427 °C. These data represent an R ratio of 0.05. Note that the [0] s experimental 
results were used to correlate the fiber damage model below 1 .04 GPa. 
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